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■ The identification of pliysical degrees of freedom is sometimes obscured in the path 

o , 

, integral formalism, and this makes it difficult to impose some constraints or to do 

' some approximations. I review a number of cases where the difficulty is overcame 

Q ' by deriving the path integral from the operator form of the partition function after 



such identification has been made. 



I. INTRODUCTION 



> 

cn 
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It is a great pleasure for me to contribute to this volume in honor of Prof. Yu. A. 

^ : 

O ■ Simonov. This is an occasion for me to remember again the origin of our friendship which 

extended over the years well beyond what it would appear from our joint papers. 
^1 The subject of my contribution is focussed on my recent interests in several problems 



which have a common feature: The identification of degrees of freedom in a path integral. 



I Indeed there are many situations where such identi ication is helpful or necessary. An old 
■ example is the thermodynamics of gauge theories but recently I met with many others. 
The first I will discuss here is how to find actions exactly equivalent to the standard ones 
but closer to the continuum at finite lattice spacinglal. This includes the definition of the 
couplings of the chemical potential, an issue of particular importance in QCD p||. The 
chemical potential is used to fix the expectation value of some charge operator. Fixing 
instead the value of the charge, namely selecting a specific charge sector is somewhat more 
difficult, but interesting pieces of information can be obtained by an expansion in a given 
charge sector of the fermion determinant in series of the number of fermions ^ . The last 
problem I will consider is how to describe the low lying excitations of fermionic systems, 



both re 
them 



ativistic and nonrelativistic, by means of effective bosons, in short how to bosonize 
. All these problems require the identification of the relevant degrees of freedom. In 
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the last case to determine the structure of the composite bosons in terms of the fermionic 
constituents, in the first case to show the equivalence of different actions. 

Now the path integral formalism is widely used because of its flexibility and the possi- 
bility of numerical applications, but the identification of degrees of freedom is not always 
easy, while it can be conveniently achieved in the operator form of the partition function. 
Therefore in all the above problems I will first identify the relevant degrees of freedom in 
the Fock space where the partition function is defined, then I will introduce the appropriate 
constraints or approximations and finally I will derive the constrained or approximated path 
integral. 

As 1 said all the subjects 1 mentioned have a common feature in the role played by 
the identification of degrees of freedom, but are otherwise very different. Therefore the 
motivations for their investigation are given separately in the relative Sections. In Section 
2 I will report the results I will use later about the standard derivation of the path integral 
from the operator form of the partition function. In Section 3 I will show how to derive 
an action different from the standard one and closer to the continuum in the nonrelativistic 
case. In Section 4 I will carry out the corresponding derivation for relativistic field theories, 
confining myself to the couplings of the chemical potential. In Section 5 I will discuss the 
case of a given charge sector. In Section 6 I will present a general method of bosonization 
valid for relativistic and nonrelativistic theories and in Section 7 my conclusions. 

II. THE STANDARD DERIVATION OF THE EUCLIDEAN PATH INTEGRAL 

FROM THE PARTITION FUNCTION 

Let me introduce some definitions. 1 denote by r the temporal lattice spacing, by Nq the 
number of temporal sites, by Xo or t the temporal component of the site position vector x, 
by T the temperature, by /i the chemical potential, by Q the (electric, baryon...) charge 
operator, and by T{xo) the transfer matrix. 

In the nonrelativistic case T is expressed in terms of the Hamiltonian 

r = exp(-T//), (1) 

and the Hamiltonian is the generator of continuous time translations. In relativistic theories 
instead only the transfer matrix is known in general, and the above equation can be used to 
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define a Hamiltonian, but only as the generator of discrete translations by the time spacing 
r. Another important difference is that the nonrelativistic interactions are generally quartic 
in the fields while the relativistic ones are quadratic. Both features contribute to make the 
derivation of the euclidean path integral different in the two cases. 

T is defined in terms of particle-antiparticle creation-annihilation operators c^, d^, c, d 
acting in a Fock space. It depends on the time coordinate Xq only through the dependence 
on it of other fields (for instance gauge fields). In fact the creation and annihilation operators 
do not depend on xq. They depend on the spatial coordinates x and on the internal quantum 
numbers (Dirac, color and flavor indices in the case of QCD), comprehensively represented 
by I, J.... 

In the transfer matrix formalism often one has to do with quantities at a given (Euclidean) 
time Xq. For this reason we adopt a summation convention over spatial coordinates and 
intrinsic indices at fixed time. So for instance, for an arbitrary matrix M, we will write 

C^M{Xo)c= Yl 4,/^x,/;y,j(Xo)Cy,J. (2) 

In this notation the charge operator Q has the form 

Q = etc - d^d. (3) 

The grand canonical partition function can be written as a time-ordered product 

Z = Tr|exp(|g)nn^o)}, (4) 

which, using the relation T^^ = tNq, and assuming the conservation of Q is conveniently 
rewritten as 

Z = Tr|n[t(xo)exp(/irQ)]|. (5) 

riQ " 

The standard way ja, l7[ to obtain the path integral form of Z is to write all the operators 
in normal order and introduce between the factors in Eq.(0) the identity 

T = J [dc^dc dd~^dd] exp{—c^c — d^d)\cd){cd\, (6) 

where the basis vectors are coherent states 

\cd) = \exp{—cc^ — dd"^)). (7) 
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The c"*", c, (i"*", d are holomorphic/Grassmann variables and satisfy periodic/antiperiodic 
boundary conditions in time for bosons/fermions respectively They have the label of the 
time slice where the identity operator is introduced. For the other indices they are subject 
to the same convention as the creation and annihilation operators. The main property of 
coherent states is that they are eigenstates of the annihilation operators 

c\cd) = c \ cd). (8) 

To get the functional form of the partition function it is only necessary to evaluate the matrix 
elements {cidi\T\cd). This, as anticipated in the Introduction, must be done in different 
ways for nonrelativistic and relativistic theories. 

The first case appears more difficult because of the quartic interactions, but since in 
general there are no ultraviolet divergencies (a detailed discussion of the departure from the 
standard form of the path integral in the presence of singular potentials can be found in ^) 
we can make without any error the approximation 

exp {-tH) r^l-rH. (9) 

Let us then consider a many-body system with the Hamiltonian written in normal form 

H{d\ a) = X! {^x^x,yfiy + flx^y'Wx.ya-ya-x} • (lO) 
x,y 

The standard expression of the euclidean path integral associated to this Hamiltonian is 

Z = y [da* da] exp(— S*) (11) 

where 

S = rY. {<+iVia, + H{al_,,, a,)} (12) 
t=i 

is the action and I denoted by 

Vt/i = -(/m-/0 (13) 
r 

the right discrete time derivative. Notice the time splitting between the fields and their 
conjugates, which implies a departure from the classical expression not only in derivative, 
but also in potential terms. Needless to say, neglecting this time splitting introduces finite 
errors jj6|. 
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In the relativistic case the transfer matrix can be written [3| 

= Tl,T., (14) 

where 

T = exp (cf Mc + d^M^d) exp (c d) (15) 

where the upper script T means transposed and I do not need to specify the matrices M, N. 
The matrix elements of T can be exactly evaluated yielding the standard form of the 
Euclidean path integral. 



III. ANTINORMAL ORDERING: A DIFFERENT ACTION IN THE 

NONRELATIVISTIC CASE 

In the standard form of the action the time is split in the fields and their conjugates. 
This is an artifact which makes the equations unnecessarily different from the continuous, 
somewhat more complicated and somewhat confusing. In gauge theories, for instance, such 
time splitting introduces a coupling of the chemical potential to the temporal gauge fields 
which has been considered of physical significance, while it can be altogether avoided in 
a different derivation of the path integral. I will first illustrate such a derivation in the 
nonrelativistic case IJ]. 

Instead of the normal order I write the Hamiltonian in antinormal order (all the annihi- 
lation operators to the left of the creation ones) 

H = ^[ho + ttx/ix.yOy + axay^^x,yaj,a],| . (16) 

where 

ho = N^{~ah^^^ + w) 
h'^y = (T/ix,y - wS^^y. (17) 

In the above equations 

w = o-t;x,x + Xl^x,z, (18) 

z 

A^^ is the number of spatial sites and a = —1 for fermions and +1 for bosons. I assume that 
the potential is sufficiently regular for w to exist, otherwise a regularization is needed or a 
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more drastic change in our procedure. Then I expand the transfer matrix and in each term I 
insert the identity between the creation and annihilation operators. For the rightmost factor 
before taking the trace I move the creation operators to the leftmost position. In this way 
the term with the time derivative remains unchanged, but in the other terms all the fields 
appear at the same time 

S' = TY,{ho + al^^Vtat + {h' - h)alat +H{a;,at)} . (19) 

t 

We can easily check on solvable models that this action gives the right results and that the 
terms arising from the rearrangement in antinormal ordering cannot be neglected. 

In the case of fermions the time splitting between the fields and their conjugates in the 
potential terms can be avoided in the same way as for bosons, but also in a simpler way. 
In fact the Grassmann fields, unlike the holomorphic variables, are independent from their 
conjugates, so that the simple transformation 

<+i - < (20) 

eliminates the time splitting everywhere, with the obvious exception of the term with the 
time derivative, which is changed into the left one. 

Because in this case we have two different derivations of the path integral, by their 
comparison we can get nontrivial identities. 



IV. ANTINORMAL ORDERING: A DIFFERENT COUPLING OF THE 
CHEMICAL POTENTIAL IN RELATIVISTIC FIELD THEORIES 

In relativistic field theories the artificial time splitting affects only the coupling of the 
chemical potential. In the Hasenfratz-Karsch-Kogut formulation j^, which is the standard 
one, for Wilson fermions such coupling takes the form 

5S = 27^5: g {[exp/i - 1] Pi,"-^ U^Ti,"-^ + [exp(-^) - 1] Pt^Tt^ul,-"^} q (21) 

X 

where q is the quark field, K is the hopping parameter, Uq the temporal link variable and 



Pt^ = ^ (1 ± 7o) 



To^^V(a:o) = /(xo±l). (22) 
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Because of gauge invariance, in the presence of the time sphtting a couphng with the tem- 
poral hnks is needed, and this led to the conclusion that a nonvanishing contribution of the 
chemical potential must necessarily involve a Polyakov loop. But Creutz showed in a toy 
model^^ that this is not true, and I will report in the sequel how to avoid such artifact 
in the full fledged QCD. 

I write the exponential of the charge in the following way 

exp(/i aoQ) = J [dc^dc dd^dd] exp{5S — c^c — d^d)\cd) {cd\ (23) 

where 

6S = {1- cosh(/iao)) (c+c + d+rf) + sinh(/iao) (c+c - d+rf). (24) 

The expression of 6S is obtained by expanding the exponential of the charge operator, 
putting all the terms in antinormal form, inserting in each monomial the unity between the 
set of annihilation and the set of creation operators and replacing them by their Grassman- 
nian eigenvalues. For the rightmost exponential of the charge before taking the trace one 
has to move the creation operators to the left of all the operators appearing under trace. 

After this, the construction of the path integral proceeds in the standard way, and we 
get the standard action with the exception of the coupling of the chemical potential where 
all the fields appear at the same time 

55* = ^ g[ (1 — cosh(/iao)) + sinh(/iao)7o]-Bq'. (25) 

X 

Here I only need to say that the matrix B does not depend on Uq. What is important is to 
notice that the quark field and its conjugate are at the same time, and then the temporal 
Wilson variable disappeared. 



V. EXPANSION OF THE FERMION DETERMINANT IN THE NUMBER OF 
FERMIONS IN A GIVEN CHARGE SECTOR 

The use of the chemical potential is a way to impose a given expectation value for some 
conserved charge. The alternative option of selecting a given sector of the charge in the path 
integral presents additional difficulties, but something can be learned by a series expansion 
of the fermion determinant 5|. To be concrete we will refer to the case of QCD, but the 
method can be applied to other cases with appropriate modifications. 
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In the absence of any condition on the baryon number the quark determinant is 

det Q = J [dqdq]exp Sq, (26) 

where Sq is the quark action and Q the quark matrix. My strategy is to write det Q as the 
trace of the transfer matrix acting in the quark Fock space, impose the restriction to a given 
baryonic sector, and then rewrite the trace as the determinant of a modified quark matrix. 
The round trip is done by mapping the Grassmann algebra generated by the quark fields 
into the Fock space following the construction of Liischer j^j . But while his paper is based 
on the mere existence of the map, in order to enforce the projection I will make use of a 
concrete realization by means of coherent states. 

The first step is then to write the unconstrained determinant as a trace in the Fock space 

detg = Trf. (27) 

The second step is to impose the restriction to a sector with baryon number rig by inserting 
in the trace the appropriate projection operator 

det Q\nB = Tr (tPub) J [dx'^ dx dy"^ dy] exp (-x+x - y+y) (x,?/|TP„^| - x, -y) (28) 

which will be expressed in terms of the determinant of a modified quark matrix. The kernel 
(x, y\'fPng \ — X, —y) has the integral form 

{x,y\'fPng\—x,—y) = J [dz~^dz dw^dw] exp (^—z^ z — w^w"^ {x,y\T\z,w){z,w\Png\—x,—y). 

(29) 

The expression of the kernel (x, ?/|T|2;, w) will not be reported here, while that of P^^ 
can easily be derived 

OO -| 

{z,w\Pr,g\-x,-y) = Y.i-ir-{{yw+r{xz+Y--+^\xx+r-+^iyyy) (30) 

Since x~^,y^ are creation operators of quarks and antiquarks respectively, we see that the 
r-th term of this series gives the gauge-invariant contribution of valence quarks plus r 
quark-antiquark pairs. 

Needless to say, for ra^ = 0, det Qlng does not reduce to the unconstrained determinant: 
Indeed also baryonic states are present in the unconstrained determinant, while they are 
absent in det Q\nB=o- In QCD at nonvanishing temperature it makes a difference whether 
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we impose or we do not the condition ub = 0. In view of the relatively low value of the critical 
temperature with respect to the nucleon mass, however, we do not expect significant effects 
from the restriction to a given baryon sector unless we go to exceedingly high temperatures. 

Let me now proceed to derive our final result. By evaluating the vacuum expectation 
values appearing in the last equation we express the kernel of the projection operator in 
terms of Grassmann variables only 

CO -| 

{z,w\P^,\-x,-y) = (31) 
To evaluate the integral of Eq. ()29j) I rewrite the above equation in exponential form 

{z,w\PnB\x,y) = E ^ g^UB+r ^"7 exp{-J,Z+X - j2W+y)\j,=,,=o- (32) 

The integrals of Eqs. ()28|) . (pi?j) are Gaussian and we get the constrained determinant in terms 
of the determinant of a modified quark matrix 

1 QriB+r Qr 

det gin. = E ^ Q^UB^r g-7 det (Q + 5Q, + 5Q,) |,,=,,=o. (33) 

The explicit form of the variations 5Q2 of the quark matrix is not important here, but 
we warn the reader that there is an error in their expression in ref.P]. 



VI. BOSONIZATION IN MANY-BODY SYSTEMS AND RELATIVISTIC FIELD 

THEORIES 

The low energy collective excitations of many-fermion systems can be described by effec- 
tive bosons. Well known examples are the Cooper pairs of Superconductivity, the bosons 
of the Interacting Boson Model of Nuclear Physics, the chiral mesons and the quark pairs 
of color superconductivity in QCD. In all these cases the effective bosons are generated by 
attractive interactions, but effective bosons can arise also in the presence of repulsive forces, 
like in the Hubbard model Some of the effective bosons are Goldstone bosons, and then 
there is a general theory which tells that they live in the coset space of the group which 
is spontaneously broken and dictates how they are related to the original fields [l^. But 
there is no general procedure to reformulate the fermionic theory in terms of the effective 
bosonic degrees of freedom, even though there are several recipes for specific cases which 
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are reviewed in [l3| . A more flexible approach is based on the Hubbard-Stratonovich trans- 
formation which linearizes the fermionic interaction by introducing bosonic auxiliary fields 
which are then promoted to physical life. The typical resulting structure is that of chiral 
theories [ij]- But in such an approach an energy scale emerges naturally, and only excita- 
tions of lower energy can be described by the auxiliary fields. Moreover in renormalizable 
relativistic field theories like QCD, the fermion Lagrangian is quadratic to start with, so that 
the Hubbard-Stratonovich transformation cannot be used. One can add quartic interactions 
as irrelevant operators, and this can help in numerical simulations, but has not led so far to 
a formulation of low energy QCD in terms of chiral mesons. 

I present a new approach p to bosonization which does not suffer from the above lim- 
itations and can be applied to theories with quartic and quadratic interactions as well. It 
is based on the evaluation of the partition function restricted to the bosonic composites 
of interest. By rewriting the partition function so obtained in functional form we get the 
euclidean action of the composite bosons from which in the nonrelativistic theories we can 
derive the Hamiltonian. In this way I derived the Interacting Boson Model from a Nuclear 
Hamiltonian. In the case of pure pairing, I reproduce the well known results for the exci- 
tations corresponding to the addition and removal of pairs of fermions, as well as for the 
seniority excitations which are inaccessible by the Hubbard-Stratonovich method. Indeed 
at least in this example this theory does not have the structure of a chiral expansion. 

For the relativistic case an investigation is in progress . 

Let me start by defining the composites in terms of the fermion operators c 



In the above equation m represents all the fermion intrinsic quantum numbers and position 
coordinates and J the quantum numbers of the composites. I assume all the structure 
matrices Bj to have one and the same dimension which I denote hj 2Q. The fermionic 
operators have canonical anticommutation relations, while for the composites 




(34) 




(35) 



It is then natural to require the normalization 




(36) 
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A convenient way to get the euclidean path integral from the trace of the transfer matrix 
is to use coherent states of composites. Therefore I introduce the operator 

j db*db{{b\b))-^\b){b\ (37) 

where 

|6) = |exp(6 -6^)). (38) 

I adopted the convention 

b-p^j:bjb^j. (39) 

J 

If the 6's where operators of elementary bosons V would be the identity in the boson Fock 
space. I would like V to be the identity in the fermion subspace of the composites. To see 
the action of V on composite operators let us first consider the case where there is only one 
composite with structure function satisfying the equation 

B^B = ^1. (40) 



Then we find 



and 



where 



m-'={^ + ^blb) (41) 



(6i|(6t)")=C„(6T, (42) 



Now we can determine the action of V on the composites 



^I(^^r)=(l-Sr7l-^) lib'T) (44) 



which shows that V behaves approximately like the identity with an error of the order of 
n/Q. It is perhaps worth while noticing that in the limit of infinite Q we recover exactly the 
expressions valid for elementary bosons, in particular 

{bi\b) = (^1 + ^blb^ ^ eMKb), n^oo. (45) 

It might appear that the treatment of states with n ~ 1] is precluded, but this is not 
true. Indeed if we are interested in states with n = n + u ior an arbitrary reference state n, 
we redefine V according to 

Pw^^Po. (46) 
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We then have 

r.\(ir)-{i-j^yii-^fwm (47) 

which shows that Vn behaves hke the identity in the neighborhood of the reference state up 
to an error of order z//(f2 — rl), namely the measure {b\b)~^ is essentially uniform. 
In the general case of many composites we have 

{h\b) = [det(l + /3i/3)]^ (48) 

where the matrix j3 is 

(3 = b-B^ (49) 

and 

{h\ibir-{bir) = ^...^exp{-Trln[l + (x.5t)(6*.5)]}U=o. (50) 

We must now make an assumption which replaces Eq. EHl namely that all the eigenvalues 
of the matrices BjBj are much smaller than Q. Then we find again that V approximates 
the identity with an error of order 

Now we are equipped to carry out the program outlined at the beginning. The first step 
is the evaluation of the partition function Zc restricted to fermionic composites. To this 
end we divide the inverse temperature in Nq intervals of spacing r 

^ = Not (51) 

and write 

Z, = tr {VTf' (52) 

where T is the transfer matrix. In the nonrelativistic case T is expressed in terms of the 
Hamiltonian 

r = exp {-tH) , (53) 

and the Hamiltonian is the generator of continuous time translations. In relativistic field 
theories instead only the transfer matrix is known in general, and the above equation can 
be used to define a Hamiltonian, but only as the generator of discrete translations by the 
time spacing r. 

At this point we must evaluate the matrix element {bi\T\b) and to do this we must 
distinguish between relativistic field theories and many-body systems. In the first case the 
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transfer matrix is a product of exponentials of quadratic forms in the fermion operators 
and the matrix element can be directly and exactly evaluated In the second case one 
must at an intermediate stage expand with respect to the time spacing r. This does not 
introduce any error because one can retain all the terms which give finite contributions in 
the limit r 0. We will report here only the nonrelativistic calculation. The most general 
Hamiltonian can be written 

H = c^hoc-J2gKl ctF^ct i cFkc, (54) 

where K represents all the necessary quantum numbers. The single-particle term includes 
the single-particle energy with matrix e, any single-particle interaction with external fields 
described by the matrix Ai and the chemical potential 

ho = e + M - jJi. (55) 

Therefore we will be able to solve the problem of fermion-boson mapping by determining the 
interaction of the composite bosons with external fields. Assuming for the potential form 
factors the normalization 

tr(Fj^Fx) =2n (56) 



and setting 



we get the euclidean action 



Tt = {\ + (3l(3t.,)-' (57) 



S = rj] j^tr[ln(l + /3;A)-lnri]-i/i + ^5^^;^tr[tr(ri/5;4)tr(r,F^A-i) 

-2 tr (r^FlF^) - ii[Ttl3;Fi, TtFKl3t-i]+] + ^tr [(5; (A_i + h A_i)] | , (58) 

where [..,..]+ is an anticommutator. This action differs from that of elementary bosons 
because 

i) the time derivative terms (contained in the first line) are non canonical 

ii) the coupling of the chemical potential (which appears in h) is also noncanonical, since 
it is not quadratic in the boson fields 

iii) the function F becomes singular when the number of bosons is of order fi, which 
reflects the Pauli principle. 
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We remind the reader that the only approximation done concerns the operator V. There- 
fore these are to be regarded as true features of compositeness. 

The bosonization of the system we considered has thus been accomphshed. In particular 
the fermionic interactions with external fields can be expressed in terms of the bosonic terms 
which involve the matrix Ai (appearing in h) and the dynamical problem of the interacting 
(composite) bosons can be solved within the path integral formalism. Part of this problem is 
the determination of the structure matrices Bj. This can be done by expressing the energies 
in terms of them and applying a variational principle which gives rise to an eigenvalue 
equation. 

The Hamiltonian can be derived by standard procedures. 

VII. CONCLUSION 

I showed that there is a number of problems which can be easily dealt with in the operator 
form of the partition function. Only afterwards it can be given the functional form of the 
path integral which is more convenient for many purposes. 

There are other examples I left over for different reasons. Among these I would like to 
mention the problem of the restriction of gauge theories to physical states, which has not 
yet found a general and satisfactory solution. 
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